Point defect concentrations in metastable Fe-C alloys 
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Point defect species and concentrations in metastable Fe-C alloys are determined using density 
functional theory and a constrained free-energy functional. Carbon interstitials dominate unless 
iron vacancies are in significant excess, whereas excess carbon causes greatly enhances vacancy 
concentration. Our predictions are amenable to experimental verification; they provide a baseline 
for rationalizing complex microstructures known in hardened and tempered steels, and by extension 
other technological materials created by or subjected to extreme environments. 
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Many industrially significant alloys are intentionally 
processed as metastable microstructures comprising a su- 
persaturation of crystal defects in various forms of aggre- 
gation How the defects are distributed and in what 
local concentrations they exist are fundamental questions 
which are not yet resolved experimentally, yet it is this 
microstructural complexity that governs the performance 
of the material. Hardened steels arc an important exam- 
ple where deformation behavior is intrinsically coupled 
to the lattice defects: in this body-centered cubic 
(bcc) Fe matrix with carbon content in excess of equilib- 
rium values for ferrite (solid solution of carbon in bcc Fe) 
and cementite (FesC), as well as high dislocation density 
and supersaturation of vacancies Q . While it is generally 
known that carbon binds strongly to crystal defects such 
as vacancies, what is lacking is a robust computational 
strategy capable of linking the fundamental physics of 
crystal defects and their energetics with quantitative, mi- 
croscopic details of the defect microstructure as it evolves 
during processing and service. 

In this work we determine the concentration of point 
defects and defect clusters in Fe-C alloys using an ap- 
proach based on first-principles calculation of the forma- 
tion energies of specific defects, and a free-energy for- 
mulation allowing either the carbon or the iron vacancy 
concentration to be out of equilibrium. We find that 
the vacancy concentration in the form of carbon-vacancy 
clusters increases dramatically in the presence of carbon 
due to strongly exothermic clustering reactions. Nev- 
ertheless, the concentrations of carbon-vacancy clusters 
remain orders of magnitude below the concentration of 
interstitial carbon. Only under the assumption of a sig- 
nificant supersaturation of vacancies will carbon- vacancy 
clusters begin to dominate the defect spectrum, in ac- 
cord with experimental observation of dominant carbon- 
vacancy clusters in irradiated bcc-Fc alloys Q ■ Ab-initio 
results on the structure and energetics of vacancies, sclf- 
interstitials and carbon interstitials in bcc-Fe have been 



reported previously |j, |5|j as well as the interaction of 
one and two C interstitials with the Fe vacancies |^. 
Our results at OK are in close agreement with cluster 
formation energies reported previously . Additionally, 
we have considered clusters with higher carbon contents 
and included di- vacancies in our considerations. 

The supersaturation of carbon in such alloys has been 
established through several experimental approaches, al- 
though the degree of supersaturation in the bcc-Fe ma- 
trix ranges from 0.3 at. % C upon 400°C annealing 
to up to 5at.% C in volumes including ferritic grain 
boundaries 0, 0|. While experiments can quantify lo- 
cal carbon concentrations (e.g., three-dimensional atom 
probe lal) or determine the carbon- vacancy binding en- 
ergy [lOj (a value of 0.85 eV was reported by positron an- 
nihilation measurements Q), it is not yet feasible to mea- 
sure the relative concentrations and specific structures of 
C-vacancy clusters in bcc-Fe during thermal processing. 
The most likely C-related defects are: (1) simple intersti- 
tial supersaturation ; (2) co-association with vacancies 
as inferred from reduced vacancy diffusion and carbon 
segregation to vacancy rich regions (3) direct obser- 
vation of elevated carbon concentration up to 8 at.% near 
dislocation cores (4) up to 6.6 at.% at grain bound- 
aries 013. 

We have performed a set of total energy calculations 
on the interaction of carbon with single and double va- 
cancies and applied the results in a statistical mechanics 
model to obtain parameter-free estimates of the concen- 
tration of point defects in various structural configura- 
tions in the bcc-Fe matrix. We employed density func- 
tional theory (DFT) [llElEl using Blochl's projector 
augmented wave method |l4L Il5l | as implemented in the 
VASP code [11 [l3 with a plane-wave cutoff of 400 eV. 
The calculations were performed in 128-atom supercells 
with the theoretical lattice constant of 2.83 A, using a k- 
mesh of 2 X 2 X 2 and a Methfessel Paxton Fermi-surface 
smearing parameter of 0.05 eV No symmetry con- 
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straints were imposed. Screening calculations and vibra- 
tional frequency calculations were carried out in 54-atom 
supercells. The geometry optimization was terminated 
with a force cutoff of 5meV/A. All calculations included 
spin polarization, starting with a ferromagnetic charge 
density. Because a wide variety of different bulk and 
defect structures are considered, we anticipate an error 
of 0.1-0.2 eV for the defect formation energies. At the 
target temperature of 160°C, where the body centered 
tetragonal phase transforms to bcc, this implies error in 
species concentration of one to two orders of magnitude. 
However, the overall conclusions regarding the dominant 
defect concentrations arc not found to be affected by this 
level of uncertainty. 

Table ^ summarizes the formation energies of the de- 
fects and defect clusters under consideration, 

-/iFe • A7lFe - ^J-C ' Anc, (1) 

where E^{T) and E°{T) denote the DFT energies of the 
supercells with and without defects present, and nx and 
Anx are the chemical potential of species X and the 
difference in atom number between the two supercells. 
Temperature dependence is taken into account through 
configurational and vibrational free energy contributions, 
the latter being calculated using the dynamical matrix 
of all atoms affected by the defect formation. To ap- 
proximate the vibrational free energy we consider each 
atom as a three-dimensional harmonic oscillator with 
force-constants derived from ab-initio total energy calcu- 
lations. The vibrational frequencies of carbon are found 
to be strongly dependent on the defect geometry, varying 
between 7 and 34 THz. Variation of the frequencies of Fe 
atoms in the vicinity of a (double-) vacancy with respect 
to the bulk value is less than 3 THz and is therefore ne- 
glected. The vibrational contributions to the formation 
energies can amount to 0.2 eV at 160°C. 

TableQlshows the most stable configuration for a given 
stoichiometry. While there are different ways to arrange a 
given number of C atoms around an Fe (double-)vacancy, 
we find that these typically have very different formation 
energies. A more extensive list of energetics of one and 
two carbon atoms near a single vacancy can be found in 
Ref. . Here we regard structures with formation ener- 
gies within 0.1 eV (our approximate error) to be degen- 
erate. Structures with energies greater than than 0.3 eV 
relative to the lowest-energy configuration arc not con- 
sidered. 

Figure^shows the stable geometries for an Fe vacancy 
surrounded by up to six carbon atoms. In the case of 
di- vacancies, we need to distinguish between two possi- 
ble orientations [100] and [111], which are energetically 
practically degenerate in the carbon-free state (compare 
Fig.lSand TableHl). Around a [111] vacancy there are two 



TABLE I: Formation energies for the different defect stoi- 
chiometries at 160° C. The chemical potentials t-tpe/c (com- 
pare Eqn. are chosen to represent thermal equilibrium of 
bulk Fe and FeaC. Xi and Xy refer to an interstitial and a 
vacancy of species X, respectively. Parentheses (...) denote 
a defect cluster. The crystallographic directions [ . . . ] refer 
to the orientation of the iron double- vacancies. 
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FIG. 1: Atomic structure of C interstitials and C-vacancy 
clusters in bcc-Fe. Fe atoms in open circles, carbon atoms 
filled circles. Left panel: octahedral and tetrahedral carbon 
interstitials; right panel: position of carbon atoms around 
a vacancy. A carbon-vacancy cluster with n carbon atoms 
contains the C atoms with labels from 1 to n. 



types of adsorption sites: those with one of the six coor- 
dinating Fe atoms missing (shaded sides in Fig. [Sjl and 
those with two missing coordinating atoms. We find that 
configurations with two missing coordinating Fe atoms 
are energetically unfavorable. Otherwise, the stable car- 
bon geometries are analogous to that of the isolated va- 
cancy. In the case of the [100] oriented di-vacancy, the 
configuration of a single carbon atom situated between 
the two vacancies is most favorable; however, this pref- 
erence does not hold for pairs and triplets of C atoms. 
From screening a wide variety of C geometries, including 
chains and tctrahcdra connecting two vacancies, no new, 
low-energy carbon configurations arc found. 

The formation of a C-vacancy cluster by combining a 
vacancy with n carbon interstitials. 
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FIG. 2: Atomic structure of C-double vacancy clusters ori- 
ented in [100] (left panel) and [111] (right panel) directions. 
Refer to Fig. for labeling. In case of the [100] vacancy, the 
position '(1)' is only occupied if there is a single C atom. 



Fey + n-C'i — y {Fey + n-C), 



(2) 



is seen to give a negative reaction enthalpy AH for n < A 
(Table QJ. The same is also predicted for two single 
vacancies/ C- vacancy clusters forming a di-vacancy/C-di- 
vacancy cluster (e.g., 2(Fey + 2C) (2Fey + 4C) [111] 
with AH — — 0.6eV). However, one should note that 
the mass-action law for a clustering reaction such as 
X + X ^ X2 relates the concentration of the cluster 
CX2 to the square of the concentration of the individual 
defects, (cx)^- Thus the defect complex is the dominat- 
ing species only if exp(— -^y) is larger than 1/cx- For 
typical defects in solids, cx is often very small (<^ 1 ppm) 
and AH is on the order of — O.leV, so this condition is 
not necessarily fulfilled. Nonetheless, enthalpy provides 
a measure of "thermal stability" of a complex once it is 
formed, as in the case of vacancy agglomerates Q. 

The equilibrium concentrations of defects Ci in a 
macroscopic crystal are given by minimizing the free en- 
ergy F{T, fipei ^J■c) sA, given temperature T and chemical 
potentials fipe and fic- 



M 

F{T,flF,,fic) = 5^C,i?f°™(T,AiFe,Mc) 

1=1 

M 

+kBT ^ [c, Inc, + (1 - c) ln(l - 



1=1 



(3) 



where E^°"^ is defined in Eqn. M denotes the number 
of different defect species and ci is the concentration of 
defect species i One obtains from Eqn.O 



Ci{T, fiFe, fJ-c) 



1 



exp 



+ 1 



(4) 



The dependence on the two chemical potentials, which 
enters through the formation energies, can be turned into 
an explicit dependence on the total carbon and/or va- 
cancy concentrations if we impose the respective con- 
straints: 



■Va/C — ^Va/C- 



(5) 



i=l 



This is useful for performing calculations under physi- 
cally motivated hypotheses. 

As discussed above, carbon and vacancy concentra- 
tions in realistic microstructures of Fe-C alloys need to 
be treated as nonequilibrium. To make the calcula- 
tions tractable we consider two alternative hypotheses 
under the condition of partial thermodynamic equilib- 
rium. That is, we have thermal equilibrium for one 
species (e.g., Fe) and non-equilibrium for the other (e.g., 
C), thus reducing the problem to a study of the effects 
of the concentration of the non-equilibrium species. In 
hypothesis (1) the carbon concentration is not in equi- 
librium with the FeaC/Fe system, whereas Fe is assumed 
to be in equilibrated with bulk Fe. This scenario is moti- 
vated by the experimental observation that carbide pre- 
cipitation is kinetically hindered and requires tempering 
at elevated temperatures 0. It is also implied that the 
Fe matrix remains substantially undistorted. In hypoth- 
esis (2) the vacancy concentration is assumed to be non- 
equilibrium, whereas carbon is assumed to equilibrate 
with the bulk Fe/FesC system. This reflects the antici- 
pated processing conditions for hardening of Fe-C alloys; 
in the course of martensitic transformation, large plas- 
tic strain leads to high intrinsic defect concentrations, 
with vacancies binding strongly to carbon interstitials. 
Due to the increased diffusion barriers associated with 
carbon-vacancy clusters Q] and the absence of a com- 
parable number of Fe interstitials (high formation ener- 
gies), it is likely that a significant vacancy concentration 
is effectively frozen in. These two hypotheses represent 
extreme scenarios; any intermediate situation can also be 
studied using the present data. 

Figure 121 shows the consequences of hypothesis (1), the 
variation of dominant defect species concentrations with 
the total carbon concentration. These results are ob- 
tained by taking the chemical potential fipe to be the 
bulk Fe value, while fic is given by the total carbon 
concentration (Eqn. EJ. We observe that the presence 
of carbon in the matrix can cause the vacancy concen- 
tration (including carbon-vacancy clusters) to increase 
by 15 orders of magnitude relative to the intrinsic va- 
cancy concentration {Fey) at a total carbon concentra- 
tion of 1 at.% Comparable increases in vacancy concen- 
tration have been reported in the presence of a different 
impurity atom, hydrogen [2llE3l. We also note that the 
concentration of carbon-free di-vacancies is negligible. 
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FIG. 3: Defect concentrations in a bcc-Fe matrix as a func- 
tion of the total carbon concentration at 160° C. Labeling as in 
Table The vertical dashed line marks thermal equilibrium 
with bulk cementite. Note that that small precipitates are ex- 
pected to have less favorable energetics in terms of atomic and 
magnetic structure, which may shift the equilibrium carbon 
concentration to higher values. 
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FIG. 4: Defect concentrations in a bcc-Fe matrix as a function 
of the total vacancy concentration at 160 °C. 



Figure 2] shows the defect concentration variation with 
total iron vacancy concentration under hypothesis (2). 
In this case [ic is fixed at the reference value, the co- 
existence of FeaC and bulk Fe. We observe that the 
number of carbon atoms associated with vacancies ex- 
ceeds the number of carbon interstitials only when the 
vacancy concentrations are larger than 10^"' at.%, which 
corresponds to a typical equilibrium vacancy concentra- 
tion near the melting point |2oj . This correlates with 
the experimental finding that C-vacancy clusters have 
been observed to be the dominant carbon-species in ir- 
radiated samples 0], which cannot be explained under 
the assumption of an equilibrated vacancy concentration 
(hypothesis (1)). 

In conclusion, our work provides the first quantitative 
assessment of the interplay of different defect species in 
Fe-C alloys. We predict defect phase diagrams as a func- 
tion of total carbon and vacancy concentrations, results 



which are amenable to experimental verification. As a 
computational framework for addressing microstructure 
complexity, our approach should be applicable to other 
advanced technological materials subject to extreme en- 
vironmental conditions. 
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